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Vortex lattice stability in the SO(5) model 
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We study the energetics of superconducting vortices in the SO(5) model for high-T c materi- 
als proposed by Zhang. We show that for a wide range of parameters normally correspond- 
ing to type II superconductivity, the free energy per unit flux T(m) of a vortex with m flux 
quanta is a decreasing function of m, provided the doping is close to its critical value. This 
implies that the Abrikosov lattice is unstable, a behaviour typical of type I superconductors. 
For dopings far from the critical value, !F(m) can become very flat, indicating a less rigid vor- 
tex lattice, which would melt at a lower temperature than expected for a BCS superconductor. 
PACS: 74.20.De, 74.20.Mn, 74.25.Ha, 74.72-h 



1. Introduction. The phase diagrams of all high- 
temperature superconductors have a rich structure, with 
two prominent features at low temperature: antiferro- 
magnetism and superconductivity. Antiferromagnetism 
(AF) is seen at low doping, while superconductivity (SC) 
is observed if the doping exceeds a critical value. 

A description of these phenomena was proposed by 
S.C. Zhang jjj, who observed that both superconductiv- 
ity and antiferromagnetism involve spontaneous symme- 
try breaking. Borrowing heavily on ideas from particle 
physics, he suggested that the symmetries involved are 
unified into a larger approximate symmetry group. He 
presented a strong case for the group SO (5), with the 
SC and AF order parameters combined, forming a fun- 
damental representation of this group. 

The parameters of the potential of the Ginzburg- 
Landau theory determine the ground state of the model; 
at high temperature, the symmetry is unbroken, while at 
low temperature, either the AF or SC order parameter 
attains an expectation value, depending on the doping. 
Because of the coupling between the AF and SC order 
parameters, exotic possibilities for solitons in the model 
can arise, as was observed already by Zhang in his origi- 
nal paper. These ideas were developed in Refs. 

In this paper, we wish to further analyze the properties 
of exotic solitons in the SO(5) model. We will consider in 
detail SC vortices [0J^), although other possibilities ||]|] 
can be analyzed similarly. We will first introduce the 
SO (5) model and review the reasons for suspecting that 
SC vortices might have AF cores. 

We will then study the free energy of vortices as a 
function of their winding number.^] Normally, in a type 



II superconductor (one for which the Ginzburg-Landau 
parameter n satisfies n > n c = l/\/2) the energy divided 
by winding number m (or energy per flux quantum) of a 
vortex is an increasing function of m. This implies that 
vortices of winding number greater than one are unsta- 
ble, which is one way of seeing that the vortex lattice is 
the preferred (lowest-energy) configuration of a supercon- 
ductor placed in an external magnetic field. For a type 
I superconductor (k < k c ), the situation is reversed: the 
energy divided by m is a decreasing function of m and 
the vortex lattice is unstable. 

However, as will be seen below, this is not necessarily 
true in the SO(5) model. Under certain circumstances, 
the vortex energy per flux of a type II superconductor can 
be a decreasing function of flux, indicating an instability 
of the vortex lattice: type I behaviour. 

The underlying reason is the possibility of an AF vor- 
tex core. When this occurs, the AF order parameter 
makes a contribution to the vortex free energy which is 
increasingly negative with increasing m. 

Two factors are involved. The first is the degree to 
which a superconductor is type II; the second is the prox- 
imity to SO (5) symmetry, which is explicitly broken away 
from "critical doping" (that which corresponds to the 
SC-AF phase boundary). These factors reinforce one an- 
other, so that a mildly type II superconductor can easily 
exhibit type I behaviour, while a strongly type II super- 
conductor requires a doping exceedingly close to critical. 

This feature of the SO(5) model gives, in principle, a 
dramatic prediction of that model. If one varies the dop- 
ing in a given superconductor, the vortex lattice should 
become less and less rigid, melting more and more eas- 



*This approach, implicit in the ground-breaking work of Bo- 
gomol'nyi M, is complementary to the usual one of studying 



the surface energy density at a boundary between normal and 
SC regions. 
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ily as critical doping is approached. Eventually, type I 
behaviour should appear. 

It must be noted that the region of parameter space 
corresponding to critical doping appears to be experi- 
mentally delicate; in particular, the appearance of in- 
homogeneities (stripe formation, phase separation) 
could mask the appearance of type I behaviour. Nonethe- 
less, reduced melting temperatures should appear away 
from this delicate region, so that an experimental signa- 
ture is still possible. Indeed, Sonier, et al. have stud- 
ied the melting of the vortex lattice in high-temperature 
superconductors and have observed melting at tempera- 
tures lower than expected in underdoped cuprates [fiof . 

2. Vortices in the SO(5) Model. According to 
the SO (5) model, the low-energy dynamics of high tem- 
perature superconductors is written in terms of a 5- 
component real field transforming as a fundamental rep- 
resentation of SO (5). The upper two components, say, of 
this real field are the real and imaginary components of 
the complex order parameter of superconductivity, while 
the lower three components are the AF order parameter. 
We will call these fields <j> = <j>x +i<j>2 and r\ — (571 , 772 , 173), 
respectively. 

The low-energy effective theory can be described in 
terms of the following free energy: 




where h = V x A is the microscopic magnetic field (hats 
will simplify notation shortly, when we go to a description 
in terms of dimensionless quantities). 

Much information (including the ground state) can be 
found by examining the potential. Including even powers 
of the fields up to fourth order, the most general potential 
is 

V{<t>,rj) = -—<)) - -T) + - (2) 

where we have written tf> = \<f>\ and 77 = \rj\. We have 
given the quadratic terms negative coefficients since this 
is what is phcnomcnologically interesting. In order for 
the potential to be bounded from below, the quartic 
terms must obey the following inequalities: 61,2 > 0, 
63 > — 

Strictly speaking, the model should be called an 
SO(3)xSO(2) model, since this is the actual symmetry 
of the model. Nonetheless, the potential is invariant un- 
der the larger group SO (5) if the two mass parameters 
are equal and if the three quartic couplings are equal. It 
will be an approximate symmetry if these couplings are 
approximately equal. In what follows, for simplicity we 



will set the three quartic couplings to the same value, 
61 = 6 2 = 63 = b. 

In order to study SC vortices, we must restrict our- 
selves to the region in parameter space that gives a SC 
ground state. This will be the case if the global mini- 
mum of the potential has a nonzero value of <fi and a zero 
value of 77. Examination of the potential shows this to 
be true if /3 = a2 2 /ai 2 < 1. Then the ground state is 
(4>, rj) = (v,0), where v = a\/\/b. It is convenient to add 
a constant a\/b to the potential, so that the free energy 
of the superconducting state in the absence of a magnetic 
field is zero. Note that /? = 1 corresponds to the SO (5) 
symmetric limit of the potential, and also to critical dop- 
ing, since neither the SC or AF state is preferred at that 
value. 

It is easy to see qualitatively why the core of a vortex 
might have an AF core (i.e., a core where 77 ^ 0). In a 
vortex (in the SO (5) model as well as in the familiar case 
of conventional superconductors), the field <f> changes in 
phase by 2tt at spatial infinity. By continuity, <f> must 
have a zero at some point, chosen to be the origin. Now 
let us look at how the field 77 fits into the situation. At 
infinity, |</5| = v and the energy is minimized for r\ = 0. 
Inside the vortex core, however, |0| — > 0. This means 
that the potential, viewed as a function of 77 with = 0, 
is minimized at 7/ 7^ 0. Were the potential energy the only 
factor, 77 would certainly develop a nonzero expectation 
value inside the core of the vortex. However potential and 
gradient energy are in competition (the gradient energy 
being minimized if 77 is zero everywhere), and the mini- 
mum energy configuration may or may not have 77 7^ in 
the core of the vortex, depending on which of these two 
competing factors dominates. The form of the potential 
suggests that as (3 is increased, there is greater likelihood 
of an AF core; this is indeed what is found numerically 
(see below, as well as Refs. J|,||]). 

As ansatz for the vortex, we use that of a conventional 
vortex (generalized to winding number m) with in ad- 
dition an ansatz for 77 (whose orientation is taken to be 
constant) which allows for the possibility of a nonzero 
core: 

4>(x) = v f( s y me (3) 

Mx) = °^^ eij s A A{s) (4 ) 
T)(x) — v n(s) (5) 

where s = r/X, A being the penetration depth, A = 
(m^^f/re* 2 *; 2 ) 1 / 2 . 

The equations of motion of the dimensionless fields 
/, n and A are (prime denotes derivative with respect 
to s): 

± (/" + \f' - (7 + ^) 2 /) + f(l-f-n 2 ) = (6) 
J_ L" + l n >^ + n tf _f_ n 2 ) = Q (7) 
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l ' + (™ +A )f 2 = Q (8) 



where in the last equation h is the dimensionless magnetic 
field, defined by h = —A' — A/s. The dimensionless 
free energy F — (2e* 2 / a\m* c 2 )F of a vortex of winding 
number m is given by 



F(m) 
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-/ 2 -/3n 2 + ^(/ 2 + n 2 ) 2 + i . (9) 



These expressions contain three parameters: the 
Ginzburg-Landau parameter k — A/£ (where the co- 
herence length is £ = (h 2 /m*ai 2 ) 1 ^ 2 ), the parameter 
(3 defined above, and the winding number of the vor- 
tex to. For high-temperature superconductors, k is usu- 
ally quite large, while (3 is determined by sample prepa- 
ration, by varying the doping. (Specifically, f3 can be 
written in terms of more physical quantities as (3 = 
1 - 8m*£(T) 2 x(Hc ~ /J 2 )/^ 2 -) P > 1 corresponds to the 
AF phase, while (3 < 1 describes the SC phase. We will 
be particularly interested in j3 < 1. 

3. Vortex energetics. For a given in, the vortex may or 
may not have an AF core, depending on the parameters 
of the model. We define /3 c (k, m), the critical value of 
(3, such that for (3 > j3 c the vortex core is AF, while for 
f3 < j3 c it is normal. Figure 1 shows (3 C as a function of 
k, for various values of m. One sees that as m increases, 
(3 C decreases. This can be understood intuitively: higher 
to corresponds to a wider vortex core, and thus greater 
impetus for 77 to attain a nonzero value in the core. 
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For a conventional superconductor, T increases or de- 
creases with to for type-II or type-I superconductors, re- 
spectively.[] The SO(5) model gives T{m) for a conven- 
tional superconductor by setting (3 — 0; then, n(s) = 
and the vortex free energy (^) is identical to that of a 
conventional superconductor. 

Figure 2 shows ^(m) for various values of (3 and k. In 
the first three plots, the upper curve ((3 — 0) represents 
a conventional superconductor: T{m) is decreasing, con- 
stant and increasing for type I, borderline I-II and type 
II superconductors, respectively. The remaining curves 
reflect the effect of an AF core in the SO(5) model. The 
fourth plot corresponds to a large value of k; (3 = is 
not displayed in order to resolve different values of (3 very 
close to 1. 

It is clear that the development of an AF core has a 
profound effect on T{m). This can be understood quali- 
tatively in the following way. As mentioned above, as m 
increases, the vortex core width increases. This is already 
true for conventional superconductors, but the effect is 
more pronounced for SO(5) superconductors when the 
core becomes AF, since in that case the free energy dif- 
ference between the AF and SC states is reduced, and 
the potential energy (which tends to reduce the core 
size) is less important. Larger core size permits a more 
spread out magnetic field, and an overall reduced energy. 
(Note that anomalously large core sizes in YBCO at low 
magnetic field have been observed though whether 
the SO (5) model can explain this has not yet been ad- 
dressed.) 

In a type I superconductor (Figure 2a) T(m) decreases 
more quickly once an AF core develops. This changes in 
a quantitative way, but not a qualitative way, the be- 
haviour of the material. 

Things are more interesting in the case of a type II su- 
perconductor (Figures 2c, 2d), where a qualitative tran- 
sition from type II behaviour to type I can be achieved. 
This occurs at approximately (3 — 0.98 and (3 — 0.9998 
for k = 7.07 and 70.7, respectively. 

Clearly for strongly type II superconductors (as is the 
case with high-temperature superconductors), (3 must be 
extremely close to 1 (doping extremely close to critical) 
for this transition to occur. Even before this point, there 
is a substantial decrease in J-(m), meaning that the en- 
ergetic savings in forming a vortex lattice (as compared 
to a large, normal region where the magnetic field pene- 



FIG. 1. /3 C as a function of k for various winding numbers. 



A very useful quantity for given values of k and (3 is the 
free energy per winding number of a vortex as a function 
of 777, J 7 (m) — F(m)/m. This quantity clearly influences 
the behaviour of a superconductor when placed in a mag- 
netic field: if T increases with m, the field will penetrate 
in vortices of winding number 1, while \IT decreases with 
to, vortices will coalesce to form large normal regions. 



^Note that we are defining the type of a superconductor 
according to the value of k rather than according to the be- 
haviour of the superconductor in a magnetic field. Since for 
conventional superconductors there is a simple relation be- 
tween the two, this distinction need not be made. In the 
SO(5) model, however, both n and j3 influence the magnetic 
behaviour of the superconductor. 
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Figure 2: J-(m) for various values of parameters k and (3. 



trates) is substantially reduced. This would be reflected 
in a less rigid, more easily melted lattice. Such behaviour 
has in fact been seen in underdoped cuprates 

We have also calculated the surface energy at a normal- 
superconducting boundary as a function of f3 and k, and 
find results consistent with the above analysis: a positive 
or negative surface energy when T(rn) is of negative or 
positive slope, respectively. This will be reported else- 
where. 

In summary, by analyzing the energy per unit flux of 
vortices as a function of winding number in the SO (5) 
model, we find that the development of an antiferromag- 
netic core has a profound effect on the behaviour of a 
superconductor in a magnetic field. This effect depends 
on the doping of the material, becoming more and more 
strong as the doping is reduced to the critical value (that 
corresponding to the AF/SC transition). More specif- 
ically, we find that the degree to which a given super- 
conductor behaves as a type II superconductor decreases 
as the doping is reduced. This can result in a less rigid 
(more easily melted) vortex lattice, and as the doping 
approaches its critical value type I behaviour results. 

Speight |l2j has recently analyzed the static intervor- 
tex force in conventional superconductivity, by treating 
the vortices as point sources. It would be interesting to 
repeat this analysis in the SO (5) model to see the effect of 
the n field on these forces, and to see if that the above be- 
haviour can be understood in terms of long-range forces 
between vortices. 

It would also be interesting to extend the work of Bo- 
gomol'nyi to the SO (5) model. This would circumvent 
much of the numerical work done in the present article. 



We have not yet succeeded in doing so, however. 
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